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ABSTRACT 



Context. Accreting black holes in galactic X-ray sources are surrounded by hot plasma. The innermost part of these 
systems is likely a corona with different temperatures for ions and electrons. In the so-called low-hard state, hot electrons 
Comptonize soft X-ray photons from the disk that partially penetrates the corona, producing emission up to ~ 150 
keV, well beyond the expectations for an optically thick disk of maximum temperature ~ 10 7 K. However, sources such 
as Cygnus X-l produce steady emission up to a few MeV, which is indicative of a non-thermal contribution to the 
spectral energy distribution. 

Aims. We study the radiative output produced by the injection of non-thermal (both electron and proton) particles in 
a magnetized corona around a black hole. 

Methods. Energy losses and maximum energies are estimated for all types of particles in a variety of models, characterized 
by different kinds of advection and relativistic proton content. Transport equations are solved for primary and secondary 
particles, and spectral energy distributions are determined and corrected by internal absorption. 

Results. We show that a local injection of non-thermal particles can account for the high energy excess observed in 
some sources, and we predict the existence of a high-energy bump at energies above 1 TeV, and typical luminosities of 
~ 10 33 erg s _1 . 

Conclusions. High-energy instruments such as the future Cherenkov Telescope Array (CTA) can be used to probe the 
relativistic particle content of the coronae around galactic black holes. 

Key words. X-rays: binaries - gamma rays: theory - radiation mechanisms: non-thermal 
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1. Introduction 

Disks around galactic black holes are formed by the ac- 
cretion of matter with angular momentum. These disks 
are opaque and emit lo cally thermal radiation of temper- 
atures up to - 10 7 K (jShakura fc Sunvaevl fl97l . These 
disks cannot be responsible for the hard X-ray emission 
extending up to a few MeV detected from Cygnus X-l 
and similar X-ray binaries and microquasars in the low- 
hard state. This finding and the presence in the spec- 
trum of features like a broad iron Ka line and a hard 
X-ray bump, led to the idea that a corona of hot plasma 
might surround the black hole and part of the disk. Two- 
te mperature models for this pla s ma w ere first suggested 
bv IShapiro. Lightman fc Eardlevl (|l976h . In these types of 
models protons are much hotter than electrons (Tj >> T e ). 
Since the pressure is dominated by the protons, the disk 
inflates, the density drops, and there is a low rate of 
Coulomb energy exchange between protons and electrons, 
allowing the existence of a two-temperature plasma in a 
self-consistent way. The Comptonization of soft photons 
results in a hard X-ray spectrum. The model, however, is 
unstable to small perturb ations in the ion temperature 
dPringldll97fit lPiranlll978r i. 
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When the plasma density is very low, the protons 
are unable to transfer energy to the el ectrons. I f mat- 
ter is advected into the black hole (jlchimarul 119771 : 
iNaravan fc Yil ll995alH) or remove d outward as a hot 
wind (jBlandford fc Begelmanl fl999). thermally stable so- 
lutions can be found. The geometry of the region with 
the hot plasma is not well-constrained, but a spherical 
region around the black hole is u s ually considered (e.g., 
lEsin. McClintock fc Naravanlll997t lEsin et al.lll998l ). The 
two-temperature plasma then forms a hot corona around 
the black hole. In the different spectral states, the cold 
disk penetrates to different distances from the black hole. 
In the very high-soft st ate it goes all the way down to 
the last stable orbit (see INaravan. Mahadevan fc Quataertl 
119981 : INaravan fc McClintock! Il998l for comprehensive re- 
views) . 

The general view of a relatively low density hot 
corona and a cold accretion disk was presented by 
iBisnovatvi-Kogan fc Blinnikovl (|1977f ). The mechanism for 
heating the corona may be mag netic reconnection of field 
loops emerging from the disk (jGaleev. Rosner fc Vaianal 
Il979j) . Violent reconnection may lead to plasma motions 
and collisions, with shock formation. A non-thermal parti- 
cle population might then arise in t he corona as the result 
of diffusive shock acceleration (e.g. JSpruitlll988l) . 

The effects of a non-thermal population of 
electrons in a hot coron a were considered by 
iKusunose fc Mineshied (|l995l) and more recently 
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bv iBelmont. Malzac fc Marcowithl f|2008h and 
IVurm fc Poutanenl ([2009T T The results of the injection of 
non-thermal protons and secondary pions and muons in a 
magnetized corona has not been comprehensively studied 
so far. The contributions from the transient particles can 
be important at high energies. The emerging emission from 
all non-thermal processes may in principle be detectable by 
future Cherenkov telescope arrays such as CTA or AGIS. 
Hence, high-energy gamma-ray astronomy can provide 
a tool to probe the non-thermal particle content of hot 
coronae around black holes. 

In this paper, we present detailed calculations of the 
radiative output of non-thermal particles in a simplified 
model of magnetized corona. The existence of the corona 
is assumed and the effect of injection of both relativistic 
protons and electrons are considered. The distributions of 
all relevant secondary particles are estimated and their ra- 
diative output computed. The coronal matter and radiative 
fields are considered as targets for the populations of rela- 
tivistic particles. The internal absorption of gamma-rays is 
calculated and the spectral energy distributions of different 
models are presented. The result is a self-consistent treat- 
ment of the non-thermal processes in the sense that once 
the medium is fixed in each model, we solve the transport 
equations for all type of particles and use the obtained par- 
ticle distributions to estimate the non-thermal radiation. 

The structure of the paper is as follows. In Sect. HJ we 
outline the basic scenario that is discussed in the paper. 
Section [3] deals with particle acceleration and losses in the 
coronal environment. The maximum energy for the different 
particles is determined. In Sect. [4[ the radiation is calcu- 
lated and the spectral energy distributions are presented, 
for different sets of parameters. Section [5] presents an ap- 
plication to Cygnus X-l. We close with a brief discussion 
and conclusions in Sects. [5] and [71 respectively. 

2. Basic scenario 

The low-hard state of accreting black holes is characterized 
by the presence of a hot corona around the compact object. 
The existence of this component is s trongly implied by the 
X-ray spectrum of Cygnus X-l (see lDovall997t lEsin et ah! 
Il998f) . Figure Q] shows a sketch of the main components of 
this system. We assume a spherical corona with a radius 
R c and an accretion disk that penetrates the corona up 
to R p < R c . For simplicity, we consider the corona to be 
homogeneous and in a steady state. 

In our model, we assume that the luminosity of the 
corona is 1 % of the Eddington luminosity of a 10 M Q 
black hole. In Table [TJ we present the values of the main 
parameters adopted and derived for the system (they are 
the typic al values of the hard X-ray luminosity of Cygnus 
X-l, see iPoutanen et all [19 97). The resulting corona has 
a luminosity L c = 1.3 x 10 37 erg s , in accordance with 
observations. 

X-ray binaries with well-formed coronae in the low-hard 
state produce relativistic jets. In these systems the power 
of the jet is understood to be related to the magnetic field, 
since the jet launchin g mechanism is likely of magnetic ori- 
gin (e.g.. lSpruitlll988[ ). For Cygnus X-l and similar systems, 
the jet kinetic power is of approximately the same luminos- 
ity as the corona. Hence, the value of the mean magnetic 
field B can be estimated by equipartition between the mag- 
netic energy density and the bolometric photon density of 



the corona L c (e.e.. iBednarek fe Giovannellill2007l and ref- 
erences therein). 
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which implies a value B — 5.7 x 10 5 G. 

It is assumed that the corona consists of a two- 
temperature plasma, with an electron temperature 
T e ~ 1 9 K and an ion te mperature Tj — 10 12 
K (e.g., iNaravan &~Yl Il995alfbh . The physical con- 
ditions in which this assumption is valid were ex- 



tensively discussed by Naravan fc Yil (|l995al fbh and 
INaravan. Mahadevan fc QuataertT i 19981 ). The corona and 
the base of the jet are thought to be a region with sim- 
i lar properties in the standard j et-disk symbiosis model 
([Malzac. Belmont fc Fabianll2009l ). We assume a spherical 
corona where thermal and non-thermal components are co- 
spatial. If the magnetic field launches the plasma into a jet, 
equipartition between magnetic and kinetic energy densi- 
ties is a reasonable assum ption that allows us to estimate 
the pl asma density (e.g., IZdziars kl 119981: INaravan fc Yil 
ll995alTbTl to be 
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where rii and n e are the ion and electron densities, respec- 
tively. For a corona mainly composed of hydrogen, this im- 
plies that rii ~ n e — 6.2 x 10 13 cm~ 3 . 

The hard X-ray emission of the corona is characterized 
by a powerlaw with an exponential cutoff at high energies, 
as ob served in several X- ray binaries in the low-hard state 
fe.g.. iRomero et aT1l2002l ) for which 



(E) = A ph E- a e- E ^ E " erg- 1 cm- 3 . 



(3) 



In ac cordance with the we ll-studied case of Cygnus X-l 
(e.g.. IPoutanen et al.l ll997h. we adopt a = 1.6 and E c = 
150 keV. The normalization constant A p h can be obtained 
from L c , 
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/ En ph {E)dE= / A ph E 1 - a e- E/E *dE. (4) 
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3. Particle acceleration and losses 

We now consider the interaction of locally injected relativis- 
tic particles with the matter, photons, and magnetic fields 
of the corona and the disk, which are taken as background 
components. 

There are three relevant processes of interaction of rela- 
tivistic electrons and muons with these fields: synchrotron 
radiation, inverse Compton (IC) scattering, and relativis- 
tic Bremsstrahlung. For protons and charged pions, there 
are also three important processes: synchrotron radiation, 
proton-proton (or pion-proton) inelastic collisions, and pho- 
tohadronic interactions. The neutral pions have a short 
mean lifetime of 8.4 x 10- 17 s, and they therefore decay 
before interacting. 

We also estimated the energy losses caused by magnetic 
Bremsstrahlung. This contribution mainly falls at interme- 
diate energies. Since this energy range is completely dom- 
inated by the thermal emission from the star and the ac- 
cretion disk, we do not show the magnetic Bremsstrahlung 
luminosity in the SEDs shown below. 
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Fig. 1. Schematic representation of the components of the system discussed in the text. In the spherical corona, the 
thermal and non-thermal components are co-spatial. Not to scale. 



Table 1. Model parameters 



Parameter 



Value 



Mbh: black hole mass [Mg] 
R c : corona radius [cm] 

R p /R c : ratio of inner disk to corona radius 

D(l + S): Covering fraction of the corona 

around thermal disk times the seed 

photon correction factor 

T e : electron temperature [K] 

Tii ion temperature [K] 

E c : X-ray spectrum cutoff [keV] 

a: X-ray spectrum power-law index 

tj: acceleration efficiency 

B c : magnetic field [G] 

rii,ne: plasma density [cm -3 ] 

A p h: normalization constant [erg 3,/5 cm -3 ] 

a: hadron-to-lepton energy ratio 

kT: disk characteristic temperature [keV] 

v. advection velocity [c] 



10« 

5.2 x I0 7(1 ' 2) 
0.08 (1) 



10 9 

I0 12 

150 

1.6 

ICT 2 

5.7 x I0 5 

6.2 x 10 13 

2.6 x I0 12 

1 - 100 

0.1 

0.1 



Typical value for Cygnus X-l in the low-hard state 
fPoutanen et al. 1997) . 

35i?G, Rg 



(2) 



GM 



The synchrotron cooling rate for a particle of mass m, 
charge e, and energy E in a region of magnetic energy den- 
sity Ub is 



t 



-l 

synchr 



4 /me V 
3 V rn ) 



In both Thomson and Klein-Ni shina regimes, the IC cooling 
rate for an electron is given by (|Blumenthal k Gould[l97fl 



e Je min Je 



YE, 

i+r 



t \ dN J 



(6) 



Here e is the energy of the incident photon, e± is the energy 
of the scattered photon, and 



dN _ 2-Krlmlc 5 n p h(e)de 



dtdei 



(7) 



where n p h(e) is the density of target photons, tq the classi- 
cal radius of the electron, and 



F(q) = 2g]n«+(l + 2 g )(l-g) + i(l-«)i^, 



r 



4eE e /{m e c 2 ) 2 , 
ei/[T(B e -ei) 



(8) 



We consider two target photon fields: the power-law 
photon field of the corona and the field from the disk. 
The latter can be represe nted by a black body with tem- 
perature kT = 0.1 keV (|Poutanen et al.lll997D . The ra- 
diation field in the corona is diluted to account for the 
solid angle subtended by the disk as seen from the corona. 
This is performed by means of a parameter D that in- 
dicates the fraction of the light emitted by the thermal 
region, that passes through the corona, and the param- 
eter S that is the ratio of intrinsic seed photon produc- 
tion in the corona to the seed photon luminosity injected 
from outside. Since we ap ply our model to Cyg nus X-l, 



(5) 

we follow the estimates of IPoutanen et all (|l997f l for this 
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source that are based on the analytic scaling approxima- 
tions for the thermal C omptonization spectrum given by 
iPietrini fc Krolikl (|1995l ). We consider the value of the dif- 
ferent parameters that are relevant to the low-hard state. 
In particular, D(l + S) = 0.08, and the ratio of the disk 
inner radius to the corona radius is taken to be ~ 0.9. For 
additi onal details, the reader is referred to iPoutanen et aD 
(Il997|). For a general picture of the C omptonization pro- 
cess, Done, Gierlihski fc Kubotal (12001 is a very useful ref- 
erence. 

For a completel y ionized plasma , the Bremsstrahlung 
cooling rate is (e.g.. lBerezinskiilll990D 



t 



Br 



4:fiiZ 2 r?,ac 



In 



2E e 
m e c 2 



(9) 



The cooling rate for inelastic collisions of protons with nu- 
clei of the corona is given by: 



VP 



pp , 



(10) 



where K pp is the total inelasticity of the process, of value 
~ 0.5. The total cr oss-section a pp can be approximated by 
(|Kelner et alJ l20061 



a inel (E p ) = (34.3 + 1.88L + 0.25L 2 



Eth 
E„ 



where 



L = In 



E n 



1 TeV 



(11) 



(12) 



77' kin 



The proton threshold kinetic energy for ir° production is 
280 MeV. 

The photomeson production takes place for photon en- 
ergies above eth ~ 145 MeV (measured in the rest-frame of 
the proton). Near the threshold, a single pion is produced 
per interaction; at higher energies, the production of mul- 
tiple pions dominates. In our model, the relevant photons 
come from the corona and, to a lesser extent, from the ac- 
cretion disk. 

The cooling rate due to photopion production for a pro- 
ton of energy E p in an isotrop ic photon field of density 
riph(e) is given by ()Steckerlll968h 



m^c 5 
~2Ej 



do 



2-,p 



2e7p 



the direct production of electrons or positrons. The cooling 
rate is also given by Eq. (|13|) , and the corresponding cross- 
section and inelasticity. The cross-section for this channel 
-also known as the Bethe-Heitler cross-section- increases 
with the energy of the photon. Both the cross-section and 
inelasticity approximati ons in the limits of low and high 
energies can be found in lBeeelman. Rudak fc Sikoral |l990) 
(see Appendix A). 

The mean lifetime of charged pions in their rest-frame 
is 7V = 2.6 x 10 -8 s, and then they decay into muons and 
neutrinos; the mean lifetime of muons is r M = 2.2 x 10 -6 
s. They decay yielding neutrinos/anti-neutrinos and elec- 
trons/positrons. In the observer rest-frame, the decay rate 



f- 1 - — 

Mec — 



tE 



(14) 



We consider two types of corona. One is an ADAF-like 
corona, where matter is advected to the black hole. This 
model wa s discussed in de tail for Cygnus X-l by iDovd 
(|1997l) and lEsin et all ()1998[ ). In this case, particles fall onto 
the compact object at a mean radial velocity v = 0.1c , 
the free-fall velocity (jBegelman. Rudak fc Sikoral [l990). 
Therefore, the convection rate is 



Rc 



(15) 



The other model considered her e is a static coron a (e.g., 
supported by magnetic fields, see lBeloborodovi[l999l) where 
the relativistic particles can be removed by diffusion. In 
the Bohm regime, the diffusion coefficient is D(E) — r g c/3, 
where r g — E/(eB) is the giroradius of the particle. The 
diffusion rate is 



L diS ~ 



2D(E) 
R 2 r ' 



(16) 



The maximum energy that a relativistic particle can at- 
tain depends on the acceleration mechanism and the differ- 
ent processes of energy loss. The acceleration rate t^ c = 
E~ 1 dE/dt for a particle of energy E in a magnetic field B, 
in a region where diffusive shock acceleration takes place, 
is given by 



de'a PJ (e')K^(e')e', (13) 



f- 1 = 

acc 



rjecB 



(17) 



where e' is the photon energy in the rest-frame of the 
proton and K m is the inelasticity of the interaction. 
lAtovan fc Dermerl (|2003l ) introduced a simplified approach 
to treat the cross-section and the inelasticity, which can be 
written as 



and 

K pl (e) 



340 /ibarn 
120 /ibarn 



200MeV < e' 
e' 



< 500MeV 
> 500MeV, 



0.2 200MeV < e' < 500MeV 



0.6 



e' > 500MeV. 



At energies below the threshold for photomeson pro- 
duction, the main channel of proton-photon interaction is 



where rj < 1 is a parameter that characterizes the efficiency 
of the acceleration. We fix 77 = 10~ 2 , which describes the 
efficient acceleration by shocks with v s ~ 0.1c in the Bohm 
regime. 

Figure [5] shows the cooling rates for different energy-loss 
processes, together with the acceleration and escape rates, 
for each type of particle considered. Under the physical con- 
ditions previously described, the main channel of energy 
loss for electrons is synchrotron radiation. Only for low- 
energy electrons IC losses are significant. For protons, both 
pp and P7 interactions are relevant. While diffusion has al- 
most no effect on particle distributions, advection plays a 
decisive role in the behavior of protons: in the model with 
advection, most protons fall onto the black hole before ra- 
diating their energy. 
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It is possible to estimate the maximum energy achieved 
by the electrons equating 
i _ i-i 



t^ 1 — t^ 1 4- / 
^acc — "-synchr ' L 



1 

IC ' 



(18) 



whereas for protons 

r 1 = r 1 +t~ 1 + T 1 

"acc v pp 1 j?7 1 "esc 



(19) 

where i csc is the timescale during which the relativistic 
particles escape from the system. This timescale is given 
by Eq. (fT5|) for models with convection, and by Eq. (TT6")) 
for models with diffusion. The maximum energies obtained 
by electrons and protons are i?max ~ 7.9 x 10 9 eV and 

Bmlx ~ 8.0 x 10 14 eV, respectively. These values are com- 
patible with the Hillas criterion, given the size of the corona. 

For pions, the main channel of energy loss is the n-f in- 
teraction, but an important fraction of these pions decay 
before cooling (those of lower energies). Muons with ener- 
gies above ~ 10 13 eV cool mostly by synchrotron radiation 
in models with a static corona. The most energetic muons 
fall into the black hole, in all models with dominant con- 
vection. 

The steady state particle distributions N(E) can be 
derived from the solution to the transport equation 
(|Ginzburg fc Svrovatskiilll964D 



_d_/dE 
dE\~dt 



loss 



N(E)) + N(E)(t- B l + t£) = Q(E) 



(20) 



where Q{E) is the injection function. The corresponding 
solution is 



N{E) 



where 



r(E,E') 



dE 




dt 


" f 

loss J E 



Q{E')e- T ( E - E ">dE', 



E' 



dE'' 



dE" 



dt 



loss 



(*csc + ^dcc)' 



(21) 



(22) 



The injection function for non-thermal protons and elec- 
trons is a powerlaw of the energy of the particles Q(E) = 
Q E~ a e~ E / Em ^, as a consequence of the diffusiv e part i- 



cle acceleration by shock waves. Following iDrurvl <| 1983 
we adopt a standard index a — 2.2. The normalization 
constant Qq can be obtained from the total power in- 
jected in relativistic protons and electrons, L rc i = L p + L e . 
This power is assumed to be a fraction of the lumi- 
nosity of the corona, L re i = q IC \L c , with q IC \ — 10~ 2 
(|Blandford fc Eichlerl fl987h . The way in which energy is 
divided between hadrons and leptons is unknown, but dif- 
ferent scenarios can be taken into account by setting L p = 
aL e . We consider models with a — 100 (proton-dominated 
scenario, as for Galactic cosmic rays) and a = 1 (equipar- 
tition between both species). 

To obtain the pion d istribution, we use t he pion injec- 
tion due to pp collisions (iFletcher et al.lll994D and jry inter- 
actions (|Atovan fc Dermerll2003l) . The injection of muons is 
calculated using the expressions from iLipari et al.l (|2007f) . 
The target fields are those of the corona in the case of mat- 
ter, and the corona-plus-disk for photons. Once the domi- 
nant energy loss process is identified, the transport equa- 
tions can be decoupled without a significant loss of accu- 
racy (say less than ~ 10%). This approximation simplifies 
considerably the calculations. A full treatment is being de- 
veloped by the authors and preliminary results have been 
used for comparison and accuracy estimates. 



4. Spectral energy distributions 

4.1. Radiative processes 

Expressions to calculat e the synchrotron spec t rum c an be 
found, for example, in iBlumenthal fc Gouldl (|1970t ). The 
power radiated by a single particle of energy E and pitch 
angle a is given by 



-fsynchr (#1 -^7 



sin ot E-y 
kmc* E-J^r 



#5/3 (23) 



where #5/3 (£) is a modified Bessel function and the char- 
acteristic energy is 



E c 



3 ehBsma ( E 
4ir mc \ mc 2 



(24) 



For a particle distribution N(E) and a volume of the emis- 
sion region V, the total luminosity is 



dEN(E)P synchl (E,E 1 ). (25) 



L-y(E^) — Ej I d 



According to IBlumenthal fc Gouldl (|1970f ). the spec- 
trum of photons scattered by an electron of energy E e = 
r y e m e c 2 in a target radiation field of density n p h(e) is 



P (F F A - 3g T c(m e c 2 ) 2 n ph (e) 
ri C (E e , £ 7 , e) — — b [q), 



(26) 



where F(q), q, and T are given by Eq. © (in this case 
61 = Ej). The allowed range of energies for the scattered 
photons is 



r 

e < £L < E e . 

- 7 - 1 + r 

The total luminosity can then be obtained from 



(27) 



L lG {E 1 )=E 2 / d 3 r I dE e N e (E e ) 



V c J E 



deP lc . (28) 



This equation takes into account the Klein-Nishina effect 
on the cross-section with energy. 

The most important products of ^7 interactions at high 
energies are pions. Neutral 7r-mesons decay with a proba- 
bility of 98.8 % into two gamma rays 



7T ^7 + 7. 



(29) 



To estimate the spectrum from the d ecay of 7T°, the 6- 
funct ional formalism can be applied (|Atovan fe Dermerl 
120031 ). The emissivity of neutral pions in this approximation 
is 

Q { ^\E^) = 5^(5^)^(5^)^0(5^) , (30) 
where uj p1 is the collision rate given by 
,^ P h(e) 



LO pl {Ep) 



de- 



2cE 



m 2 p e > 
2E 2 J Hh. 



pip 



de'a$(e')e', (31) 
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Fig. 2. Radiative losses in a corona characterized by the parameters of Table [TJ 



and n^o is the mean number of neutral pions created per 
collisi on. For more details, we refer to lAtovan &; Dermerl 

pooah . 

Taking into account that each 7r° decays into two pho- 
tons, the photon emissivity is given by 

<7 7 (£ 7 ) = 2 J Q%?\E v )6(E y - 0.5E„)dE„ 

= 20iVp(10£ 7 )cjp 7 ^(10£^)n 7r o(10£ 7 ) ■ (32) 

It is possible to estimate the injection function of pions pro- 
duce d by pp interactions also usin g the ^-functional formal- 
ism {Ah aronian fc Atovanl I2000T ) . In this approximation, 
the emissivity of 7r° is 

Q% P) {E«) = cn, J S(E n — n 7T Eki n )a pp (Ep)Np(E p )dEp 

= <T P p I rripC H I N p I m p c H I (.66) 



For proton en ergies in the GeV-TeV range, k s» 0.17 
(|Gaisserlll990t ) and the total cross-section a pp can be ap- 
proximated by Eq. ([XT]) for E^ > 1 GeV, and a pp « for 
E kin < 1 GeV. 

Once the pion injection Q 7r o(i? 7r ) is known, the gamma- 
ray emissivity q~ f (E~ f ) (erg -1 cm -3 s -1 ) can be obtained 
from 

^= 2 Ljw=^ dE " <34) 

where E min = E 1 + m^c^/AEy. 
4.2. Secondary pair injection 

Two channels for secondary pair production were taken into 
account: pairs produced in photon-photon interactions and 
electron and positron production as a result of muon decay. 

The energy spectrum of pairs produced in photon- 
photon interactions was studied, for example, by 
lAharonian et al.l (|l983l ). Under the conditions e -C m e c 2 < 
.Ely, the pair injection Q e (^e) (in units of erg _1 s~ x cm -3 ) 
can be approximated by the expression 
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Qe(E e 



3 COT 

32 m e c 2 



de^dui 



n 7 (e 7 ) riph(w) 



4e 2 

7 



4 TeU 7 -7e) 



7e (e 7 - 7 e ) 



In 



2(2e 7 ^ - l)e 7 _ ( ^ 



5e 7 o; + 



7 W / 7 e 2 (e 7 - 7e) 2 



(35) 



where 7 e = E e /m e c 2 is the Lorentz factor of the electron, 
e 7 = E~ ( /m e c 2 , and cj = e/m e c 2 are the dimensionless pho- 
ton energies. 

The second production mechanism of secondary elec- 
trons/positrons is the decay of charged pions, given by 



7T ± -> + ^(^)> 



-> e ± + v^U/j) + v e (y e ). 



(36) 
(37) 



This process has been extensively studied by several au- 
thors. Given a muon distribution N^j^), it is possible to 
obtain the pair injec tion Q e (7e) as a result of this process 
using ()Ramatvlli974l ) 



s(7e 



dry. 



where 7^ max = 104, 



it = 7e7e ± V7e ~ 1V7? - ^ 



(38) 



(39) 



and the electron distribution in the rest-frame of muon is 



Pile, 



v /max 3 



(40) 



4.3. Spectral energy distributions 

Figure [3] shows all contributions to the total luminosity for 
different advection regimes and for different values of the 
parameter a. On the one hand, it can be seen that the 
luminosities produced by primary electrons are higher in 
models with a = 1. On the other hand, luminosities pro- 
duced by hadrons and muons are higher in models with a 
static corona and diffusion of the relativistic particles. This 
is because in models with advection an important fraction 
of protons and pions are swallowed by the black hole, while 
with diffusion these particles are able to lose their energy 
before falling onto the compact object or escaping from the 
system. In models with a = 100, the non-thermal emission 
at > 1 MeV is dominated by synchrotron and IC radia- 
tion of secondary pairs. At very high energies, the main con- 
tributions to the spectrum are due to photo-meson produc- 
tion in all models. We note that below ~ 150 keV the source 
will be totally dominated by thermal Comptonization (not 
shown in the figures for clarity). 

From Fig. [3j we conclude that there are two param- 
eters that determine the relevant radiative processes: the 
hadronic content in the plasma and the advection velocity. 



If the hadronic content is high, then a large number of sec- 
ondary particles are expected to increase the emission at 
high energies. This is precisely what happens in our model 
for a high value of the parameter a. However, the advec- 
tion also has an important role, because in a corona with 
convection a significant part of the proton content will be 
engulfed by the black hole reducing the emission. The over- 
all SED predicted by a particular model is then the result of 
the specific balance between the two main free parameters. 

There are other physical quantities that are important 
to our model, such as the magnetic field and the background 
thermal luminosity, but the values of these parameters are 
restricted by observations of sources, such as Cygnus X-l. 

4.4. Absorption 

Gamma-rays created in astrophysical sources can be ab- 
sorbed by different mechanisms. The absorption can be 
quantified by the absorption coefficient or opacity r. If the 
original gamma-ray luminosity is Ly(Ej), the attenuated 
luminosity L~ I {E 1 ) after the photon travels a distance I is 



L 7 (£ 7 ) = L°JE^)e- T 



(l,E-,) 



(41) 



The opacity depends on the absorption process. In the 
model presented here, the main mechanism of absorption 
is photon-photon pair production. This interaction takes 
place between gamma-rays and thermal X-rays photons 
from the corona. This process is possible only above a kine- 
matic energy threshold given by 



E^e > (m e c 2 ) 2 , 



(42) 



in a frontal collision. 

For a photon created at a distance r from the center 
of the corona, the opacity cause d by photon-photon p air 
production can be obtained from (jRomero fc Vilall2008h 



'77 



n p h(e, r')<7 77 (e, E~f)dr'de, 



(43) 



where e is the energy of the X-ray photons and n p h their 
densi ty. The total cross-sec tion is given by the expression 
(e.g., lAharonian et al.|[l985| ) 



/3 2 ) (3-/3 4 )ln 



1-/3 



where 



(me£T^ 1/2 
E 7 e 



+ 2/3(/3 2 - 2) ,(44) 



(45) 



Figure |4] shows the difference between the produced 
luminosity and the attenuated luminosity. For 10 MeV 
< Ej < 1 TeV, almost all photons are absorbed. 

Although the gamma-ray distribution depends on the 
electron/positron distributions and viceversa, we treat each 
population independently. This assumption is justified be- 
cause the strong magnetic fi eld is understood to prevent 
the d evelopment of cascades (|Pellizza. Orellana Sz Romero! 
I2010h . 

The main cooling channel for high-energy elec- 
trons/positrons is synchrotron radiation, since the IC in- 
teraction occurs mainly in the Klein-Nishina regime. This 
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Fig. 3. Spectral energy distributions obtained for different sets of parameters. Internal absorption not included. 



can be clearly seen from the first panel of Fig. [5J The syn- 
chrotron cooling then halts t he development of high-energy 
casca des in the corona (see iBosch-Ramon fc Khangulvanl 
120091 ). Only at low energies (£ 7 < 100 MeV) may a third 
generation of pairs be produced (jVievro et al.ll20l6T) . 



5. Application to Cygnus X-l 

We applied the model discussed in the previous sections to 
the well-known binary system Cygnus X-l. Cygnus X-l is 
a very bright X-ray binary consisting of a compact object 
of ~ 10.1 Mp) and a companion 09.7 lab star of ~ 17.8 M Q 
(iHerrero et al.lll995[). at an estimated distance of ~ 2 kpc 
(e.g.. iGierlihski et al.lll997l and references therein). The X- 
ray emission alternates between soft and hard states. The 
spectrum in both states can be approximately represented 
as the sum of a bl ack body and a pow erlaw with expo- 
nential cutoff fe.g.. IPoutanen et al.l Il997h . During the soft 
state, the black body component is dominant and the pow- 
erlaw is steep, with a photon spectral index ~ 2.8 (e.g., 
iFrontera et al.l [20011 ). During the hard state, more energy 



is in the powerlaw compon ent, which is even hard ened, with 
photon indices ~ 1.6 (e.g., Gicrliriski et al.l[T997[ ). 

iMcConnell et al.l ([2000) reported a high-energy tail 
in the low X-ray state of the source, extending from 
50 kcV to - 5 MeV. The data at MeV energies, col- 
lected with the COMPTEL instrument of the Compton 
Gamma-Ray Observatory, can be described as a power- 
law with a photon spectral index of 3.2. Observations with 
the INTEGRAL satellite have confirme d the existence of 
a sup ra-thermal tail in the spectrum (jCadolle Bel et aD 
2006 ). So-called hybri d ther mal/ non-the r mal m odels from 
IPoutanen fc Svenssonl (|!996t ) and lCoppil (1 1999ft have been 
used to fit the observed spectrum ([McConnell et aD 1*2000: 
iCadolle Bel et al.ll2006[ ). These models consider a hybrid 
pair plasma with a non- thermal component. In particular, 
using the EQPAIR code. lMcConnell et al.l (|200Q[ ) concluded 
that either the magnetic field in Cyg X-l is substantially 
below equipartition (at least to within an order of magni- 
tude) or the observed photon tail has a different origin than 
that related to locally accelerated electrons. 

We applied our lepto-hadronic model to Cygnus X-l. 
In Fig. [SJ we show the predictions of the models discussed 
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Fig. 4. Produced luminosity and attenuated luminosity for different sets of parameters. 



in the previous sections (equipartition and proton domi- 
nance) to Cygnus X-l. All our models assume equiparti- 
tion magnetic fields. As expected, the emission in the MeV 
range is dominated by products of hadronic interactions 
and secondary pairs. The best fits are for a model with 
a = 100 and little or null advection, with absorption play- 
ing a major role in shaping the spectrum. At high-energies 
Ey > 1 TeV, a bump produced mainly by photo-meson 
production appears. It might be easily detectable by the 
future Cherenkov Telescope Array, if Cygnus X-l is within 
its declination range. In Fig. [5] we indicate the sensitivity of 
different instruments, including MAGIC, Fermi, and CTA. 



MAGIC detecte d a high-energy flare from Cygnus X-l 
(|Albert et al.ll2007l) . which, ho wever, was likely produce d 
in the jet of the object (e.g., iBosch- Ramon et al"1l2008h . 
Similar consider ations apply to the fl are detected by the 
AGILE satellite (|Sabatini et al.ll2010j) . In our model, even 
if large magnetic reconnection events were to modify the 
non-thermal population on short timescales, the GeV emis- 
sion would be totally suppressed by photon annihilation in 
the thermal bath of the corona. Gamma-ray flaring events 



at GeV energies, then, cannot arise from a strongly mag- 
netized corona. 



6. Discussion 

Proton-proton interactions in a purel y thermal non- 
magne tize d corona were considered by Eilek fc Kafatod 
(| 19831 ) and lMahadevan. Naravan fc Kroiild (|l997tl . The lat- 
ter authors also considered the effects of a non-thermal 
population of protons, but without taking into account 
pair pro duction in the soft photon field, which was consid- 
ered by Jourdain fc Rogues (I1994T) and independently by 
iBhattacharvva et al.l (|2003l l2006t ). The effects of the mag- 
netic field were not considered by any of these authors. 

The role of the magnetic field in the cooling of the dif- 
ferent types of relativistic particles in the corona is very 
important and cannot be ignored, as it can be clearly seen 
in Fig. [5] Electrons are cooled almost immediately by syn- 
chrotron losses in fields of ~ 10 5 G, truncating the devel- 
opment of IC cascades. Photo-meson and photo-pair pro- 
duction in models with large values of the parameter a can 
produce significant high-energy gamma-ray emission and 
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Fig. 5. Application of the model to Cygnus X-1. Observational data from COMPTEL ()McConnell et al.l[2000l ). We 
include the sensitivity of MAGIC, Fermi, and CTA. 



secondary pair injection. The high-energy emission may be 
detectable in the future from different sources. Since it is 
of hadronic origin, detections or upper limits can be used 
to place constraints on the number of relativistic protons 
in the corona. 

Emission in the range 100 MeV to 1 TeV is not expected 
because of absorption in the soft photon field. All emis- 
sion detected in this ran ge should be produced in the je t 
(e.g., iRomero et al~ll2002l 120031; iBosch-Ramon et al.|[2008h . 
However, if a sudden injection of relativistic protons occurs, 
for instance as a consequence of major reconnection events, 
a neutrino burst may be produced. This p ossibility will 
be ex plored in detail in a forthcoming paper (jVievro et al.l 
2010). 
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Appendix A: Bethe-Heitler process and relevant 
approximations 

The cross-section for proton-photon pair production is also 
known as the Bethe-Heitler cross-section. This function in- 
creases monotonically with photon energy. In the limits 
of low an d high energy, the analytical app roximations are 
given by (Bc gelman. Rudak fc Sikoralll990D 



7. Conclusions 

We have developed a model to predict the radiative out- 
put of a two-temperature magnetized corona with a mixed 
population of non-thermal electrons and protons. All ra- 
diative processes have been calculated in this scenario. The 
contribution of secondary pairs and transient particles such 
as muons and charged pions has been shown to be impor- 
tant. The complexity of the spectral energy distributions 
calculated reflects the variety of the physical processes in- 
volved. We have applied our model to the case of Cygnus 
X-1 and obtained a good agreement with the observed soft 
gamma-ray tail observed by COMPTEL and INTEGRAL. 
We have also predicted a high-energy bump at E > 1 TeV 
that might be detectable with future giant telescope arrays. 
Additional research into the inhomogeneous coronae is in 
progress. 



r(«) 



(x') fa 1.2 x 10~ 27 



/(e)' 
th 



th 



cm 



for 2 > x' < 4, and 



(*') 



5.8 x 10 



3.11n(2x') -8.1 + 
(2.71n(2x') -ln 2 (2<) 



+ 0.71n 3 (2x') + 0.5 



cm 



(A.l) 



(A.2) 



for x' > 4. The variable x' is such that e' — x'm e c 2 is 
the photon energy measured in the proton rest-frame and 

x'fa = 2 is the threshold energy. The inelasticity of this 
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process has a maximum at x'l h , and decreases monotoni- 
cally with photon energy. The low and high energy approx- 
imations are 



, m, 



K^(x') = 4 v 



m„ 



+ 0.00781n(a;' - 1) 
for x' < 1000, and 



1 + 0.4 In (x' - 1) + O.lln (x' - if + 

3" 



(A.3) 



m v 



-8.8 + 5.6 Ins' - 1.61n:z;' 2 + OJlnz' 3 ' 
3.11n2x'-8.1 



|(A.4) 



for x' > 1000. 

The resulting pairs are then considered as a secondary 
source of leptons when estimating the spectral energy dis- 
tribution. This contribution is, in general, negligible com- 
pared to other sources of secondary pairs in the context 
discussed in this paper. 



References 

Aharonian, F. A., Atoyan, A. M., Nagapetyan, A. M., 1983, 

Astrophysics, 19, 187 
Aharonian, F. A., Kririllov-Ugrimov, V. G., Vardanian, V. V., 1985, 

Ap&SS, 115, 201 
Aharonian, F. A., Atoyan, A. M., 2000, A&A, 362, 937 
Albert J., et al. 2007, ApJ, 665, L51 
Atoyan, A. M., Dermer, C. D., 2003, ApJ, 586, 79 
Bednarek, W., Giovannelli, F., 2007, A&A, 464, 437 
Begelman, M., Rudak, B., Sikora, M., 1990, ApJ, 362, 38 
Belmont, R., Malzac, J., Marcowith, A., 2008, A&A, 491, 617 
Bcloborodov, A. M., 1999, in J. Poutanen & R. Svensson (eds), High 

Energy Processes in Accreting Black Holes, ASP Conference Series 

161, p.295 

Bcrczinskii, V. S., 1990, Astrophysics of Cosmic Rays, North-Holland, 

Amsterdam 
Bhattacharyya, S., et al. 2003, ApJ, 595, 317 

Bhattacharyya, S., Bhatt, N., Misra, R., 2006, MNRAS, 371, 245 
Bisnovatyi-Kogan, G. S., Blinnikov, S. I., 1977, A&A, 59, 111 
Blandford, R. D., Begelman, M. C, 1999, MNRAS, 303, LI 
Blandford, R. D., Eichler, D., 1987, Phys. Rep., 154, 1 
Blumenthal, G. M., Gould, R. J., 970, Rev. Mod. Phys., 42, 237 
Bosch-Ramon, V., Khangulyan, D., 2009, Mod. Phys., 18, 347 
Bosch-Ramon, V., Khangulyan, D., Aharonian, F. A., 2008, A&A, 
489, L21 

Cadolle Bel, M., et al. 2006, A&A, 446, 591 

Coppi, P. S., 1999, in J. Poutanen & R. Svensson (eds), High Energy 
Processes in Accreting Black Holes, ASP Conference Series 161, 
p.375 

Done, C, Gierlinski, M., Kubota, A., 2007, A&A, 15, 1 
Dove, J. B., 1997, ApJ, 487, 759 

Drury, L., 1983, Reports on Progress in Physics, 46, 973 

Eilek, J. A., Kafatos, M., 1983, ApJ, 271, 804 

Esin, A. A., McClintock, J. E., Narayan, R., 1997, ApJ, 489, 865 

Esin, A. A., et al. 1998, ApJ, 505, 854 

Fletcher, R. S., et al. 1994, Phys. Rev. D, 50, 5710 

Frontera, F., et al. 2001, ApJ, 546, 1027 

Gaisser, T. K., 1990, Cosmic Rays and Particle Physics, Cambridge 

University Press, Cambridge 
Galeev, A. A., Rosner, R., Vaiana, G. S., 1979, ApJ, 229, 318 
Gierlinski, M., Zdziarski, A. A., Done, C, et al. 1997, MNRAS, 288, 

958 

Ginzburg, V. L., Syrovatskii S. I., 1964, The Origin of Cosmic Rays, 

Macmillan, New York 
Herrero, A., et al. 1995, A&A, 297, 556 
Ichimaru, S., 1977, ApJ, 214, 840 



Jourdain, E., Roques J. P., 1994, ApJ, 426, L11-L14 
Kelner, S. R., Aharonian, F. A., Bugayov, V. V., 2006, Phys. Rev. D, 
74, 033018 

Kusunose, M., Mineshige S., 1995, ApJ, 440, 100 

Lipari, P., Lusignoli, M, Meloni, D., 2007, Phys. Rev. D, 75, 123005 

Mahadevan, R., Narayan, R., Krolik, J., 1997, ApJ, 486, 268 

Malzac, J., Belmont, R., Fabian, A. C ., 2009, MNRAS, 400, 1512 

McConnell, M. L., et al. 2000, ApJ, 543, 928 

Narayan, R., Yi, I., 1995a, ApJ, 444, 231 

Narayan, R., Yi, I., 1995b, ApJ, 452, 710 

Narayan, R., Mahadevan, R., Quatacrt, E., 1998, in M.A. 
Abramowicz, G. Bjornsson, and J.E. Pringle (eds), Theory of Black 
Hole Accretion Disks, Cambridge University Press, p. 148 

Narayan, R., McClintock, J. E., 1998, New Astronomy Reviews, 51, 
733 

Pellizza, L., Orellana, M., Romero, G. E., 2010, IJMP D, in press 
Pietrini, P., Krolik, J. H., 1995, ApJ, 447, 526 
Piran, T., 1978, ApJ, 221, 652 

Poutanen, J., Krolyk, J. H., Ryde, F., 1997, MNRAS, 292, L21 
Poutanen, J., Svensson, R., 1996, ApJ, 470, 249 
Pringle, J. E., 1976, MNRAS, vol. 177, 65 

Ramaty, R., 1974, in F. B. McDonald and C. E. Fichtel (eds), 
High Energy Particles and Quanta in Astrophysics, MIT Press, 
Cambridge (Mass.), p. 122 

Romero, G. E., Kaufman Bernado, M. M., Mirabel, I. F., 2002, A&A, 
393, L61 

Romero, G. E., et al. 2003, A&A, 410, LI 

Romero, G. E., Vila, G. S., 2008, A&A, 485, 623 

Sabatini, S., et al. 2010, ApJ, 712, L10 

Shakura, N. I., Sunyaev, R. A., 1973, A&A, 24, 337 

Shapiro, S. L., Lightman, A. P., Eardlcy, D. M., 1976, ApJ, 204, 187 

Spruit, H. C, 1988, A&A, 194, 319 

Stecker, F. W., 1968, Phys. Rev. Lett., 21, 1016 

Vieyro, F. L., Romero, G. E., Vila G. S., 2010 (in preparation) 

Vieyro, F. L., Romero, G. E., 2010 (in preparation) 

Vurm, I., Poutanen, J., 2009, ApJ, 698, 293 

Zdziarski, A. A., 1998, MNRAS, 296, L51 



